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^ ■ Abstract 

■ The convergence to the equihbrium of the solution of the quantum Kac model for 

Bose-Einstein identical particles is studied in this paper. Using the relative entropy 



^ I method and a detailed analysis of the entropy production, the exponential decay 

■ rate is obtained under suitable assumptions. The theoretical results are further 

_ , illustrated by numerical simulations. 

< 

: 1 Introduction 

In this paper we study the equation governing the time evolution of a gas composed of 
Bose-Einstein identical particles. Let f{t,v) be the velocity distribution function at time 
t > with the velocity f G M. According to quantum physics, the presence of a particle 
in the velocity range dv increases the probability that a particle will enter that range: 
' the presence of f{v) dv particles per unit volume increases this probability in the radio 

cn ■ 



> 



1 + f{v). Following Chapman and Cowling this fundamental assumption yields the 
^ . so-called Boltzmann-Bose-Einstein equation, that is the quantum Boltzmann equation for 

Q I Bose-Einstein particles. This equation has been extensively studied in physical literatures 
I and numerical simulations. However, there are not many rigorous mathematical results. 
We mention here for spatial homogeneous isotropic case, a theory of weak solutions devel- 
oped by Lu in [T71 (TS] , and another class of locally defined in time classical solutions by 
X : Escobedo et al. in [HI [131 E] • See [21 US] for more reviews of currently available math- 
i ematical results. On the other hand, AUemand and Toscani in |1] derived the following 
nonlinear Fokker-Planck equation (Kac model) 

dtf = Af{t)d,J + Bf{t)d,{vf{l + /)) (1) 

with 



Af{t) = I v'f{l + f)dv, Bj{t) = I fdv. 
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This model is obtained as the grazing coUision hmit of one- dimensional Boltzmann equa- 
tion for Bose-Einstein particles in the spirit of Kac caricature of a Maxwell gas with a 
singular kernel [l]. However, the existence of good solutions for this integro partial differ- 
ential equation is still unknown, and is currently under investigation. We will furthermore 
specialize to this Kac model, and study the convergence of the solution for the Kac model 
to the Rose distribution by using entropy method. 

The rigorous study of the convergence to equilibrium is by now classical in kinetic 
theory. For example, using the classical logarithmic-Sobolev inequality of Gross [15], 
and the Csiszar-KuUback-Pinsker inequality |6l HE], the convergence to the equilibrium 
with exponential decay rate can be derived by the relative entropy method for linear 
Fokker-Plank equation. For the nonlinear Fokker-Plank-Landau equation, the trend to 
equilibrium has been obtained by Desvillettes and Villani in [8]. Toscani and Villani in 
[22] studied the convergence to the equilibrium for the Boltzmann equation. Except for 
the spacial homogeneous kinetic models, Desvillettes and Villani [10\ studied the trend 
to equilibrium for the spacial inhomogeneous linear Fokker-Planck equation. For more 
about the trend to equilibrium for classical kinetic equations, we refer to [HI |9l [7] and 
references therein. 

In [1], Carrillo, Rosado and Salvarani studied a 1-D quantum Fokker-planck equation 

dtf = d,J + d,{vf{l + f). 

Note that the factor 1 + / comes from the quantum effects. It is easy to see that, the 
mass, / f{t,v)dv, is conserved along the time evolution. By using the relative entropy 
method, they proved that the solutions converge to the Bose equilibrium with exponential 
decay rate. The above model, a simplified model of (II]) with Af{t) and Bf{t) replaced by 
constant 1, does not conserve the kinetic energy. 

For the Kac grazing limit model in comparing with the model studied in there 
is an additional conservation law: conservation of kinetic energy, i.e J |f p/df. However, 
the nonlinearity in the Kac model is stronger. For later use, let m and e be the mass and 
the energy defined by the initial data fo{v), 



= y /o(^^) dv, ^ = y ^ /o(^^) dv, 
by supposing that /o > 0, /o G -^^^(M), and f^/o G L^(]R). The entropy H{f) is defined as 
H{f){t) = / l{f){tiv)dv, with 7(x) = a; logo; — (1 + a;) log(l + x). 



We further remark that the entropy used in [1] is the sum of the entropy defined above 
and the kinetic energy which is conserved for the Kac grazing limit model ([1]). 

We shall work with smooth enough solutions: we show in the next section that one 
can get a priori weighted LP' bounds on solutions, and similar estimates hold also true for 
higher derivatives. 

The main result in this paper is stated as the following Theorem. 

Theorem 1. Let f{t,v) be the solution of the Kac grazing limit model ([T]) with initial 
data fo which is positive and satisfies 

fo{v) dv = m, / v'^foiv) dv = e 
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for some positive constants m, e. Assume ||/o||^2 < m/2 and m?/e he suitably small. 
Then there exists a positive constant a depending on Jq, such that 



\\f{t)-U\L 



^<C{U)e-^'/\ 



with constant C{fo) > 0. Here foo is the Bose distribution with mass m and energy e 
defined by 

In comparing the results obtained in [3] for quantum Fokker-Planck equation, some 
additional assumptions on the initial data are needed in Theorem [H In fact, a generalized 
logarithmic-Sobolev inequality developed in [5] for nonlinear diffusion equation was used 
directly in [1] to control the entropy production from below by the relative entropy. In 
their proof, an auxiliary nonlinear diffusion equation which has the same entropy and the 
equilibrium was introduced. While for the Kac model ([1]), it is impossible to introduce 
such auxiliary equation with the same relative entropy or the equilibrium state, and a 
compatible entropy production term. Without using the generalized logarithmic-Sobolev 
inequality, we follow some ideas used in |3], we get the decay rate of entropy production, 
then the convergence of the solution of the Kac grazing limit model to its equilibrium. 
The constraints on the initial data stated in Theorem [1] will be used to get the decay rate 
of the entropy production. We believe that these constrains are only needed to simplify 
the proof. While in the numerical simulation part, we don't take into account these 
constraints. 

The rest of the paper is organised as following. In Section [2], we will give some 
preliminary estimates. Then the entropy and the entropy equality will be introduced 
in Section |3l The Bose distribution /oo will be given from the equivalent form of the 
entropy production. Then based on a detailed study on the entropy production, we get 
the exponential decay rate. Finally the theoretical results will be illustrated by numerical 
simulations in Section HI 

2 Preliminaries 

In this section we will show some a priori weighted estimates on the solution, together 
with some control on a specific quantity Af{t). These estimates will be used in next 
section in order to get the decay rate of the entropy production. 

Before starting our estimates, we note first that it is not difficult to assert the positivity 
of the solution / as in |4j. Let us repeat their arguments quickly here: let be the 
Friedrich mollifier, and define the smoothed sign and absolute functions 




Multiplying Kac equation by (pe{f), and integrating it over M, we get 




Mi) U'sif)\dJ\'dv-Bjit) vfil + fmf)djdv. (2) 
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Note that 
and 

- I vf<P'^{f)dJdv = - I vdM'Mf)-fMf))dv- I vd^l^l^ e^P'eiOdC^ dv. 

The first term on the right hand side of ((21) is non positive since 0^ > 0. And the sec- 
ond term on the right hand side of ([2]) vanishes as e — ?• from Lebesgure's dominated 
convergence theorem. Then we have by letting e — )■ on both side of ([2]), 



II/WI-/WIIli<|||/o|-/oII 



If the initial data /o G is non negative a.e on M, then the solution / (if exists) belongs 
to and is always non negative a.e on M. 

In conclusion, we have shown that if /, the smooth solution of Kac's model with 
initial data /o G i^^(M), is sufficiently decaying, there holds that the norm of / is 
non- increasing for t > 0. Furthermore if /o is non-negative a.e. in M, the solution f{t,v) 
is also non-negative in M for any t > 0. 

2.1 Weighted estimates of the solution 

In this paragraph, we are going to show the following uniform estimates in time. 
Proposition 1. The norm of the solution f verifies 



\\f(t)\\i. < max <; ||/o||i„ — + -— !> , fort>0. 

Further, assume that m^/e is sufficient small and ||/o||^2 < m/2, then there exists a 
positive constant a depending on J'q, m and e, such that 

nf{t)\\h-m<-a, fort>0. (3) 

Similarly, let < v >= (1 + |f p)"*^/^, we have the weighted estimate 

\\<v> /(t)||i2 < max{|| <v> /o||i2,i;(m,e)}, 

for t > 0, with constant Y^{m, e) depending on m and e. 

Proof. Firstly we multiply the equation by / and then integrate the resulting equality 
with respect to v to get 

Ift^imih + Mt)\\djmh = m(^l\\fmh + ^ii/wiiia) • (4) 

Then we use the interpolation inequality and the Gagliardo-Nirenberg inequality (cf. 
(3.27) in [20]) to estimate the norm of / as 

ii/wilia < \\fmL4fmu < Aii/wiu.ii/(t)iu.|i5./wiu^ 

< 4djmh + ^\\fmi^, 



where A denotes the constant arising in the Gaghardo-Nirenberg inequahty. Recall that 
Af{t) = e + 1 1 "w/ 1 1^2 • Choose £^ = ^ and finally we have 



d / A^m'^ 

-j^\\fmh+emmi,<(m+^ 



h+e\\djmh< rn+-—]\\f{t)\\l,. (5) 



The Nash inequality (cf. (6) in [21j) in one dimensional case reads 



where Cat is a numerical constant. Using the Nash inequality in ([5]) gives 
^ll/Wlli. <-C2Ke)||/(t)||i.+Ci(m,e)||/(t)||i„ 

with Ci(m, e) = m + A^m'^/(9e), C2{m,e) = Cj^e/m^. The above differential inequality 
can be solved explicitly in a standard way. For simplicity let us feel free to omit the 
dependence of the constants on the mass m and the energy e and denote X{t) = ||/(t)|||2, 
then we have 

X'it) < -C2X%t) + CiX{t), 

which can be reduced to 



Integrating the above inequality over [0,t] gives the upper bound of as 



X\t) < ^ < max X(0)^ \ . 

X(0)2 ^ Ci I ^ 

Inserting the expressions of Ci and C2 in the above inequality, we get the uniform bound 
for the norm as stated in Proposition [H 

Next, it is easy to check that when m^/e is small enough such that ^ < m/2, we get 
the existence of the positive constant a stated in Proposition [1] thus the inequality ([3]) 
holds. 

The L?' estimation of < > / is similar. Let g =< v > f with < v >= (1 + v'^Y^'^. 
We multiply the Kac model hy < v > g, then integrate it with respect to v over R to 
get 



1 d ,2 



Af{t) J <v> gd,Jdv + Bj{t) J <v>gd,ivfil + f))dv 
2 , < V >^ 3 ./ < f 



Since 1/ < v >< 1, \v\/ < v >< 1, we have 



Similarly we use the interpolation inequality and the Gagliardo-Nirenberg inequality to 
estimate the norm as 



umis < M\9mL49mL4dMt)\\L^ 



We take e — 3e/(2m) and use Nash inequality to get 

As < f; >< 1 + v'^/2, we have \\git)\\L^ < m + f . Note that Af(t) < e + \\g(t)\\l^. Then 
we have 

^^\\9ml^<Ci{m,e)\\g{t)\\h + 2\m^^^ 

with 

Ci(m,e) = 2e + m + — ^ ^' 



9e 
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Let — ||5'(i)|||,2 and the differential inequality can be written as 

Y\t) < G{Y{t)) = Ci{m, e)Y{t) + 2Y{tf - C^im, e)Y^{t). 
As G{Y^ has a unique positive zero point 

1 + a/i + C'i(m,e)C'2(m,e) 
y,(m,e) = ^- ^ > 0, 

and G'(F) is positive over ]0, Y^\ and negative on +oo[. Then we get the global existence 
of which will take values between the initial value i^(0) and the equilibrium point 
y*. In conclusion, we have 

||5/(t)||i2 < max{||^o||i2,y*}, 

with gro =< > /o- 

Remcirk 1. 1. Note that we have the following interpolation inequality 



27/2 

which follows classically by optimizing w.r.t. R>0 the following inequality 

f fdv+ f fdv<\\f\\L2V2R+^. 

J\v\<R J\v\>R -n- 



m 



Therefore an LF' control on f implies a control of ■ 

2. We have only shown weighted estimations of f but similar estimates on higher 

derivatives also hold true. It is important to note that smoothness is not required 
for estimating convergence to equilibrium. 
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2.2 Estimate for the quantity Af{t) 

In this paragraph we study the time derivative of the quantity Af{t) = e + ||u/|||2. 
Proposition 2. There holds 

A'Jt) 

<2|i/(t)ih 



A'(t) 



Proof. Let g{t,v) = vf(t,v), which verifies 

dtg = Af{t)d,,g - 2Af{t)dJ + Bj{t)d,{vg + - Bf{t)g{l + /). 

Then multiply the equality by g and integrate it with respect to v. Finally we get the L 
equality as 



2 



= Af{t)\\f{t)\\l,-B,{f) J f{t,v)g\t,v)dv. 

As Af{t) = e + \\g{t)\\l2, we have 

A'fit) + Af{t){Bf{t) - 2\\f{t)\\l) < Bf{t)e < Bj{t) ■ Aj{t). 
Using the definition of Af{t) we complete the proof of Proposition [2l 

3 Relative entropy method and decay to equilibrium 

In this section, we will prove Theorem [T] by the relative entropy method. Firstly we will 
introduce the entropy, the entropy production and the equilibrium to the Kac model ([T]). 
Secondly we will show the decay rate of the entropy production. Finally the decay rate 
of the solution for the Kac model to the equilibrium can be derived. 



3.1 Entropy and equilibrium 

Let 7(/) = / log / - (1 + /) log(l + /). Note that 7'(/) = log j^j- The entropy H{f) 

H{f) = J l{f)dv 

verifies the entropy equahty 

^^(/) = - / (^/W/ffT) + Bf{t)vdJ^ := -D{f). (6) 

The entropy production D{f) can be written in some other forms. For example since 
^/(^) = Ivf = ^ Iv ^(/) ^^"^ written as 

D{f) = J (^Aj{t)^^^^+2Bj{t)vdJ + B{tr'^ dv 



1 



Afit) 



f{l + f)\Af{t)d,j'if) + Bf{t)vfdv. (7) 



From ([ZD we get that D{f){t) > 0. 

Furthermore, we can use the expression of Af{t) into the entropy production D{f) 
and write it in a symmetric form as 

= fil + f)m + f*)\vd.,l'if*)-v*d,l'{f)\'dvdv., 

with = /(f*). From the equahty 

the equihbrium foo is 

fooiv) = / ■ , (8) 

exp(Aif^ - A2) - 1 

where the constants Ai > and A2 will be determined by the initial data. Note that the 
equihbrium f^o defined above is the so-called Bose distribution function. 



Remark 2. We can show that 



A, = (9) 

J 00 



In fact, we have 

Bj^ = ^ I : , d^; 



Ai j - 



Af = I —. 5 r—- df = — (9\, / — ; 5 —dv. 



Using the expressions of Bf, we have 



A. = -dx.Bf^. = / -dv = —Bf 



Note that ([9]) can also be obtained from the entropy production in the form 

In conclusion, we have the following lemma 
Lemma 1. The equilibrium foo minimizes 

H{f) : f{v) is positive, I f{v) dv = m, I v'^f{v) dv = 



with m and e fixed. As 7 is convex, this minimizer function is unique. Moreover, given 
any solution f{t,v) to the Kac model ([1]) with initial data fo of mass m and energy e, we 
have 

H{foo)<H{f){t)<H{fo), t>0, 

and 

hm H{f){t) = H{foo). 

r— f 00 

Before ending this paragraph, we introduce the relative entropy H{f\f 00) as 



H{f\foo) = H{f) - HUoo) = j Mf) - 7(/oo) - 7'(/oo)(/ - foo)] dv, 
where we used the conservations of mass and energy for the last equality. 
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3.2 Decay rate of the entropy production and the relative en- 
tropy 

To get the decay rate of the entropy production, we shall study the time derivative of 
D{f). To simplify notations we denote ^ = Af{t)d^'y'{f) + Bf{t)v. Hence the Kac equation 
and the entropy production D{f) can be written as 

dj = d4f{i + m], = -J^ I /(I + /)e 

Then we have 



:= I + 11 + III. (10) 

Next we will calculate these three integrals. Firstly the integral I can be written as 

A' it) 
Af{t) 

Then the second integral II can be calculated as 

We denote by 0(/) = f^ + ^P + f. Then using the expression of ^ = Af(t)d^'y'{f) + Bf{t)v, 
we can rewrite II as 

II = -Ml [ d. + 2 / ( ^'^f^ - 2^11111] \.d f P 

Aj{t) J "f"^^^^ ^^ + '7 v/(i + /) Pil + fy J'^ 

+ ^jjl^^^,idvfdv[f{l + m]. (11) 

Finally, from the conservation of mass, Bf{t) = m, we get 

A' (t) f r A' it) /"I 

III = y j f{i+mdiiu) = 2-^z^(/)-2 y y^^i5.[/(i+/)e]r. 

In summarize, we get the derivative of the entropy dissipation 
where {■ ■ ■ } denotes some positive terms. As 0(/) > /(I + /), then we get 



We use first Proposition [2] and then Proposition [T] to get 

^'^^^^ i?K^)<2||/(t)||i.-5/(t)<-«. 



Af{t) 
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In the last inequality we used the smallness assumptions as in Proposition [TJ Immediately 
we derive the following decay rate about the entropy production 

D{f){t) < D{fo)e-^\ 

We use the decay rate of the entropy production D{f) in the entropy equality to 

get 

^i^(/|/oo) > -/^(/o)e-"*. 
Then integrating the above inequality over ]ti,t2[ gives 

p—at2 p—ati 

HiflUih) - HiflUih) > D{fo) . 

a 

Let t2 — )■ +c>o, and as hm^^oo H{f\foo)(t) = 0, finally we get the decay rate of the relative 
entropy 

H{f\U{t) < DifoY- 



a 



3.3 Decay rate of the norm and the proof of Theorem [T] 

Next, we show the decay rate of the solution to the equilibrium. Observe that there 
exists a function y{t, v) which takes values between f{t, v) and fooiv) such that the relative 
entropy can be written as 



H{f\U = H{f)-H{U = /(^(/)-7(/^)-y(/^)(/-/^))d^; 

l"{y{t,v)){f-Wdv, 



Remark that we used the property of mass and energy conservations and the Taylor 
formula in the last two equalities. 

As using the Cauchy Schwartz inequality, we have 

11/(0 < 



< J /oo(i + /oo) J i"{y{t,v)){f-uy 

„-at 

< C HU\U<C DU,) . 

a 

Hence using the mass conservation we get the following desired result 
\\f{t) - /ooIUhm) = 2||/(t) - /oo||lh/</^) < C(/o)e-" 
Hence the proof to Theorem [1] is completed. 
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4 Convergence towards equilibria: numerical simula- 
tions 



Since we have shown the solution goes exponentially fast to the Bose equilibrium distri- 
bution, we will do some numerical simulations in this section, to show the equilibrium 
distributions for different initial states, and the exponential decay of the entropies. 
We recall first the Bose distribution 



fo 



1 



exp(Aif 2 - A2) 



The numerical simulations are carried out by different initial conditions. The first 
example shows if the initial data is concentrated near the center, it will evolve to Bose 
distribution, with entropy decaying exponentially to some final state. 
Example 1. Consider initial data 

0.1 



g(i'-7r/2)2+0.1 _ 

The equilibrium distribution and evolution of entropy are shown in Figure [H 



(12) 



Entropy Evolution 





200 400 600 800 1000 1200 1400 1600 1800 2000 



Figure 1: Initial data fo in ( 1T2|) . 

This property is true without smallness assumption on the initial data. In Example 2, 
we take 10 times the initial data as in Example 1 and observe also the exponential decay 
of the entropy, with different time scale used in the simulation. 
Example 2. Consider initial data 

•^0 ^ g(t;-7r/2)2+0.1 _ X' ^^'^^ 

The equilibrium distribution and evolution of entropy are shown in Figure [2l 

As we know the Bose distribution behaves like Gaussian when |f | is big. The next 
example shows the evolution of a Gaussian to Bose distribution. 
Example 3. Consider initial data 

fo = 5*e-^'/'. (14) 
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- Initial data 

- Equilibriun 



Entropy Evolution 





Figure 2: Initial data /o in ( IT3l) . 




- Initial data -9-9'^ 

- Equilibrium 



Entropy Evolution 




Figure 3: Initial data /o in ( IT4l) . 



The equilibrium distribution and evolution of entropy are shown in Figure [31 

The comparison shows the Bose distribution is more singular near |v| = , but bahaves 
like a Gaussian for |f | big. 
Example 4. Consider initial data 



fo = 8*[ 



(15) 



The equilibrium distribution and evolution of entropy are shown in Figure HI 

This example shows the evolution of summation of two Gaussians. We will show a 
more general case in next example. 
Example 5. Consider initial data 



/o 



2 + -^ for ve [-^,0], 

5 2, ,571, 
for ve[0,—], 

for others. 



(16) 
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10 I — I 1 1 1 1 1 1 — Entropy Evolution 




Figure 4: Initial data /o in ( ITSll . 
The equilibrium distribution and evolution of entropy are shown in Figure [51 



7 1 — I 1 1 1 1 — Entropy Evolution 




-15 -10 -5 5 10 15 t 



Figure 5: Initial data /o in ( !T6|) . 

Note that this initial data is not in space LlogL, since /o = in some interval thus the 
entropy is oo at first several time steps. This more general case also shows the exponential 
decay of entropy. 

All the above numerical results showed the quick convergence towards to equilibria, 
especially, we can see the exponential evolution of entropy. The numerical results further 
elaborate our main result stated in Theorem [TJ 
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